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ABSTRACT 

Using HSTAVFPC2 data, we have performed a comparative study of the Blue Straggler Star (BSS) popu- 
lations in the central regions of the globular clusters AM 1, Eridanus, Palomar 3, and Palomar 4. Located at 
distances R cc > 50 kpc from the Galactic Centre, these are (together with Palomar 14 and NGC 2419) the most 
distant clusters in the Halo. We determine their colour-magnitude diagrams and centres of gravity. The four 
clusters turn out to have similar ages (10.5-1 1 Gyr), significantly smaller than those of the inner-Halo globulars, 
and similar metallicities. By exploiting wide field ground based data, we build the most extended radial density 
profiles from resolved star counts ever published for these systems. These are well reproduced by isotropic 
King models of relatively low concentration. BSSs appear to be significantly more centrally segregated than 
red giants in all globular clusters, in agreement with the estimated core and half-mass relaxation times which 
are smaller than the cluster ages. Assuming that this is a signature of mass segregation, we conclude that AM 
1 and Eridanus are slightly dynamically more evolved than Pal 3 and Pal 4. 

Subject headings: globular clusters: general — globular clusters: individual (AM 1, Eridanus, Pal 3, Pal 4) 



1. INTRODUCTION 

In the colour-magnitude diagram (CMD) of a globular clus- 
ter (GC), Blue Straggler Stars (BSSs) define a sub-population 
located along the main-sequence (MS) in a position brighter 
and bluer than the current MS-turnoff (TO). For this reason 
these objects are thought to be hydrogen burning stars more 
massive than a typical TO star. Two physical mechanisms are 
thought to be respo nsible for their formation: mass-transfer 
in a b inary system (lMcCreal[T964l) . and direct stellar colli- 
sions (Hills & Day 1976). Since collisions are more frequent 
in regions of higher density, the relative importance of the dif- 
ferent formation channels could depend on the e nvironment 
(e.g., iFusi Pecci etail Il992t iDavies et all 120041) : BSSs in 
loose GCs might preferentially arise from mass-tr ansfer activ- 
ity in primordial binarie s (hereafter MT-BSSs; iFerraro et al.l 
120061: iLeigh et al.11201 lal) . while BSSs located in high density 
environments mi ght mainly form fro m stellar collisions (here- 
after COL-BSSs: lFerraro et al.ll2004l). 

Discovered for the first time bv lSandagel (119531) in the exter- 
nal regions of the GC M3, it was with the advent of the Hub- 
ble Space Telescope (HST) and the 8-meter class telescopes 
that it became possible to search for BSSs in dense clus- 
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ter cores (e.g. iParesce et al.l 119911 IFerraro & Parescel 119931 
IClark et aLj }2004). BSSs ha ve been found in any properly im- 
aged GC dPiotto et alj2004]) and they are studied also in dwar f 
spheroidal galaxies dMapelli et al] [2007: Mo nelli et aTJl2012l) . 

Since the BSS formation mechanisms seems to be tightly 
connected with the cluster internal dynamics, these stars are 
commonly recognized as ideal test particles to investigate 
the impact of dynamics o n stellar evolution in different envi- 
ronmental conditions (e.g. | Bailvnlll 995: Sills & Bailvnll 19991 
Moret tiet al.ll2008t IFerraro & Lanzonill2009l) . An interesting 
example comes from the dis covery of two distinct sequences 
of BSSs in the core of M30 (IFerraro et al.ll2009l) . The authors 
argued that each of the two sequences is populated by BSSs 
originated by one of the two formation channels, both trig- 
gered by th e collapse of the clus ter core a few Gyr ago. On the 
other side, Knigge et al. (2009) suggest that most BSSs, even 
thos e found in cluster cores, come from binary systems (see 
also lLeigh et al.1201 lbl) . Nevertheless the parent binaries may 
themselves have been affected by dynamical encounters. 

Systematic studies of the BSS populations in GCs have 
shown that in most cases their radial distribution is bimodal, 
with a high peak in the cluster centre, a trough at intermedi- 
ate radii and a rising br anch in the external regions (see e.g. 
IDalessandro et al. 20 091 and references therein). Dynamical 



simulations (e.g. iSigurdsson et al.lll994t iMapelli etalJ.120041 
l2006tlLanzoni et al.ll2007allbl ) suggest that the observed shape 
of the BSS radial distribution depends on the relative con- 
tribution of COL-BSSs and the efficiency of the dynamical 
friction, that progressively segregates objects more massive 
than the average (as BSSs and their progenitors) towards the 
cluster cen tre. Very interestin g exceptions are the cases of uj 
Centa uri (IFerraro et al.ll2006l). NGC 2 419 (IDalessandro et al.l 
120081) and Pal 14 dBeccarietal.ll2011l) , where the BSS radial 
distribution is completely flat. This fact indicates that in these 
clusters the dynamical friction was not effective yet in seg- 
regating BSSs toward the centre of the potential well. Inter- 
estingly enough NGC 2419 and Pal 14 are among the most 
remote GCs in the Galaxy, lying at a distance from the Galac- 
tic centre Rcc > 50 kpc. 
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FIG. 1.— (V/V-I) CMDs of the GCs AM 1, Eridanus, Pal 3, Pal 4. All 
the sampled stars are shown with black dots, while the location of the mean 
ridge line is shown with grey points. We fitted the CMD with isochrones 
from lDotter"eTaTI H2008bt grey solid lines). The location of the MS-TO is 
also indicated. 

We decided to extend the investigation to other four Galac- 
tic GCs (namely AM 1, Eridanus, Pal 3, Pal 4) located in 
the extreme outer-halo. In Section [2] we describe the data-set 
and reduction procedure. The CMD and the age estimate for 
each cluster are discussed in Section |2~T1 The determination 
of the astrometric solut ions and the cluster centres of gravity 
is presented in Section 12.21 The radial density profiles from 
resolved star counts and the estimates of the cluster struc- 
tural parameters and characteristic time-scales are discussed 
in Section[3] Section|4]is devoted to the radial distribution of 
the BSS, while a summary is presented in Section[5] 

2. TFIE DATA SET 



IStetson et ail (119991 hereafter S99) published deep opti- 
cal CMDs obtained with the Wide Field Planetary Camera 2 
(WFPC2) on board the HST for Eridanus, Pal 3 and Pal 4. 
The core region of Eridanus was roughly centred on the WF3 
chip of the WFPC2 mosaic, while Pal 3 and Pal 4 were cen- 
tredon the PC chip. We refer to S99 (see also Harris et al. 
119971) for a detailed description of data quality an d reduc- 
tion p rocedure. In brief, a standard DAOPHOTII (IStetsonl 
1987) Point Sprea d Function (PS F) fitting procedure, includ- 
ing ALLFRAME ( Stetson] [l 9941) . was adopted to obtain in- 
strumental magnitudes and colours for all measurable stars 
in the WFPC2 fields. Here we adopted the S99 photomet- 
ric catalogues of Pal 3, Pal 4 and Eridanus listing the F555W 
andF814W magnitudes of 2968, 4122 and 2340 stars, respec- 
tively. 

AM 1 was observed with the WFPC2 in the F555W and 
F814W filters during Cycle 6 (GO-651 2; PI Hessert. The 
observation log is described in Table 1 of Dotter et al. (2008a, 
hereafter D08). We retrieved the AM 1 images from the STScI 
archive and performed independent PSF fitting photometry. 
The PSF was modelled on each image using a number (be- 
tween 40 and 80) of isolated and well sampled stars, adopting 



FIG. 2. — The (V,V-I) CMDs of AM 1 (crosses), Eridanus (open squares), 
Pal 3 (open triangles), and Pal 4 (open circles) superimposed one to the other 
and aligned at the TO colour and magnitude, (V-I)to and Vjo, respectively. 

the PSF/DAOPHOTII routine. A first PSF fitting was then 
performed on every image using DAOPHOTII/ ALLSTAR, 
while ALLFRAME was used to obtain a fine measure of the 
star magnitudes. For each star all the magnitudes were nor- 
malized to a reference frame and averaged together, and the 
photometric error was derived as the r.m.s. of the repeated 
measurements. 

For homogeneity with the S99 re duction procedure, we 
used the colour terms in table 7 from iHoltzman et al.l (11995b 
to convert the instrumental (F555W,F814W) magnitudes to 
the standard Johnson V and Kron-Cousins I (hereafter V and 
I, respectively). 

2.1. Colour-magnitude diagrams and age 

The CMDs obtained for AM 1, Eridanus, Pal 3 and Pal 4 
are shown in Figure Q] (see also figures 1-3 in S99). All the 
detected stars are shown as black dots, while the location of 
the mean ridge line obtained with a second order polynomial 
interpolation with a 2.5 sigma clipping rejection criterion is 
displayed with grey points. The high quality of the WFPC2 
images allows us to properly sample the clusters' stars from 
the Tip of the Red Giant Branch (RGB) down to ~ 3 mag 
below the MS-TO. 

The solid line in Figure Q] sh ows the location of the best-fit 
theoretical isochrone from D otter et al.l (l2008bl) over-plotted 
on the CMDs of each cluster. We used the distance modu- 
lus and reddening from iHarrisI (1 1 9961 2010 edition ), while we 
adopte d m etallicity [Fe/Hl=-1.5 8 and -1.41 from lKochet alj 
(2009) and lKoch& Cdtel bOlOh for Pal 3 and Pal 4, respec- 
tively, [Fe/H]=-1.42 from S99 for Eridanus, and [Fe/H]=-1.5 
from D08 for AM 1 . With these parameters we obtain an ab- 
solute age of 10.5 Gyr for Eridanus and Pal 4, and 1 1 Gyr for 
Pal 3 and AM 1, confirming that these clusters are younger 
than the inner-halo GCs (S99 and D08). 

The location of the TO, defined as the hottest point along 
the MS in the theoretical isochrone of each cluster, is also 
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TABLE 1 

Center of gravity of the target clusters. 



Cluster 


RA 


Deco 


AM 1 


03 h 55 m 02.5 s 


-49°36'53.2" 


Eridanus 


04 h 24 m 44.7 s 


-21°11'13.9" 


Pal 3 


10 h 05 ra 31.6 s 


00°04'21.7" 


Pal 4 


ll h 29 ra 16.5 s 


28°58'22.4" 



marked in Figure [TJ This was used to compute the shifts in 
magnitude and colour that, after taking into account also the 
different distance moduli and reddenings, are necessary to co- 
add the four CMDs. The result is impressive (see Figure [2): 
the RGB morphologies of the four clusters are in full agree- 
ment and consistent with a small metallicity difference among 
their stellar populations. Indeed, the combined CMD mimics 
a highly homogeneous single-metallicity population. These 
results are in agreement with the study of S99 and D08, who 
find that these clusters are younger than M 3 by 1.5-2 Gyr, and 
they are homogeneous in terms of both age and metallicity. 

2.2. Astrometry and centre of gravity 

The coordinates of the stars identified in each cluster have 
been converted from the "local" to the absolute astrometric 
system by using a well tested p rocedure adopted by our group 
for more than 10 years (see e.g. lLanzoni et al.l2010l and refer- 
ences therein). Since only a few (if any) primary astrometric 
standards can be found in the central region of GCs, we usu- 
ally complement HST photometry with ground based wide- 
field imaging. 

B- and R-band images of AM 1 where taken in December 
2000 with the Wide Field Imager (WF£3) at the MPE/ESO 
2.2m telescope (66.B-0454(A); PI: Gallager). Pal 3, Pal 4 and 
Eridanus were observed in the g and r bands with MegaCam, 
the ~ 1° x 1° FoV CCD camera at the Canada-France-Hawaii 
Telescope (CFHT). We retrieved the bias-subtracted and flat- 
field corrected images from the CFHT Science Data Archive 
(Observing run ID: 2010AC06 for Pal 3 and Pal 4; 2009BC02 
for Eridanus; PI Cote). We analysed the data following the 
procedure described in Section [2] We analysed only the few 
chips (from 2 to 4, depending on the location of the cluster 
centre in the FoV of each data-set) which allow us to sample 
the entire cluster radial extent. Notice that the tidal radius o f 
the four clusters ranges between - 2' and - 5' dHarrisll 19961) . 
Unfortunately the ground based observations reach only the 
TO level and were not deep enough to properly sample the 
BSS region with an appropriate signal-to-noise ratio. For this 
reason we used this data set only to search for astrometric so- 
lutions and to supplement the HST data for the construction of 
the cluster projected density profiles from resolved star counts 
(see Sect. [3]). 

We used more than one hundred primary astrometric stan- 
dard stars from the GSC2.3 catalogue in the cluster vicinity to 
derive an astrometric solution and obtain the absolute equato- 
rial (RA and Dec) positions of the stars sampled in the ground 
based catalogues. Cross-correlations of the catalogues and 
astrometric solutions were calculated with CataXcorr, a code 
developed and maintained by Paolo Montegriff o at the LNAF- 
Bologna Astronomical O bservatory (see e.g. lLanzoni et al.l 
l2010tlBellazzini et alJ201 lb . Finally, the positions of the stars 
in the WFPC2 data were transformed into the same celestial 

10 The WFI is an imager composed of a mosaic of 8 CCD with a pixel 
scale of 0"238/pix and a total FoV of ~ 34' X 33'. 



coordinate system by adopting ~ 50 stars in common with 
the ground based catalogues. The r.m.s. scatter of the final 
solution was ~ 0.3" in both RA and Dec for the four clusters. 

The centre of gravity (C grav ) of the four cluste rs is estimated 
as the barycenter of the resolved stars (see e.g. lLanzoni et alj 
2010). In doing this an iterative procedure was adopted. The 
method estimates the centre by averaging the RA and Dec po- 
sitions of all the stars contained within a circular area of a 
given radius and proceeds until convergence is reached. We 
performed this procedure using three different limiting radii 
and, at any given radius, using only stars with magnitudes 
brighter than Vto, Vto + 0.5 mag and Vto+ 1 niag. The final 
value of Cgrav is calculated as the average of the nine mea- 
sures (see Table [TJi. The uncertainties both in RA and Dec 
are smaller than 1 arcsec for Eridanus and AM 1, while they 
increase to ~ 2" for Pal 3 and Pal 4 because of the extremely 
low stellar density even in the core region of these clusters. 
Within the errors the new ce ntres are in g eneral good agree- 
ment with the ones listed bv [Hanisl (fT996b . 

3. RADIAL DENSITY PROFILES 

We used the photometric catalogues and the estimated val- 
ues of Cgrav t0 calculate the projected radial density profile 
from resolved star counts for each cluster. We divided the en- 
tire sample into concentric annuli, each split into four sub sec- 
tors (quadrants). The number of stars within each subsector 
was counted, and the star density was obtained by dividing 
these values by the corresponding subsector areas. The stellar 
density in each annulus was then obtained as the average of 
the subsector densities. 

The star counts in the clusters' central regions were per- 
formed using the stars with a magnitude V < {Vjo + 1 ) from 
the WFPC2 high-resolution imaging catalogues, while the 
stars sampled down to the TO were taken from the ground 
based data set for the external regions. Notice that the stars 
adopted to compute the density profile cover a very limited 
range of mass in both the HST and the ground based datasets, 
so that a single mass-model can be used to properly describe 
the observed profile (see next Section). The final density pro- 
files of the program clusters are shown in Figure [3] where the 
star counts from ground based catalogues (open diamonds) 
have been normalised to the WFPC2 ones (solid circles) us- 
ing at least two circular area in common between the two data 
sets. These are the most accurate and extended radial density 
profiles ever published for these GCs. 

3.1. King models and structural parameters 

The derived density profiles allowed us to derive the clus- 
ters' str^turalparameters through fitting isotropic King mod- 
els dKing||196 6). The background density was estimated from 
the outer data points. The best fit has been determined by % 2 
minimization. 

As shown in Figure [3] the observed profiles are very well 
reproduced by isotropic King models characterised by the 
quoted parameters, with Wo, c, r c and r, indicating the central 
dimensionless potential, the concentration, the core and tidal 
radii, respe ctively (see also Table |2j we adopted the distances 
quoted by Hamslll996)) . These parameters have been used to 
estimate the core relaxation time (t rc ) and the half-mass relax- 
ati on time (t r h) of the four clusters from equation (10) and (11) 
of iDiorgovskil dl993l). We assumed the average stellar mass 
(m*) = 0.42Mrj~ iiSollima et all 1201 ll) while the total cluster 
masses M c \ have bee n derived from the observed total lumi- 
nosity (Harris 1996) and the mass-to-light ratios quoted by 
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FIG. 3 . — Observed surface density profiles in units of number of stars per square arcsecond. Star counts from the WFPC2 data are shown as solid circles, while 
those from ground based catalogues, normalised to the former, are shown as open diamonds. The solid line represents the King model that best fits the observed 
density profile over the entire cluster extent. Its characteristic parameters are labeled. The dashed horizontal line shows the adopted stellar density background, 
while the dotted-dashed vertical line indicates the cluster tidal radius (r,). The lower panel shows the residuals between the observations and the fitted profile. 



TABLE 2 

Age and structural parameters of the four GCs. Ages were estimated through the best-fit of theoretical isochrones. The un certainties on the 

AGE ARE OF THE ORDER OF 0.2-0.3 GYR. CLUSTER TOTAL MASSES HAVE BEEN COM PUTED FROM THE OBSERVED TOTAL LUMINOSITIES QUOTED By[HARRIS I 1996]) AND 
THE MASS-TO-LIGHT RATIOS ESTIMATED BY lMCLAUGHLIN & VAN DER MARLL I 2005 ) (AN AVERAGE VALUE OF 1 .9 HAS BEEN ADOPTED FOR THE MASS-TO-LIGHT RATIO OF 
ERIDANUS, WHICH IS NOT INCLUDED IN THIS LATTER WORK). THE DIMENSIONLESS POTENTIAL, CONCENTRATION, CORE AND TIDAL RADII AND CENTRAL MASS DENSITY 
(Wo, c, r c , r, AND log(po), RESPECTIVELY) ARE FROM THE BEST-FIT KING MODELS. THE CORE RELAXATION TI ME l0g(« re ) AND HALF-MASS RELAXATION TIME log(f rft ) ARE 

COMPUTED FOLLOWING EQ. (10) AND (11) OF DlORGOVSKI I 1993), RESPECTIVELY. 



Cluster 


Age 


log(M d ) 


W 


c 


r c 


n 


log(po) 


log(f rc ) 


log(^) 




[Gyr] 


[Mq] 






[pc] 


[pc] 


[M Q /pc 3 ] 


[yr] 


[yr] 


AMI 


11.0 


4.10 


5.60 


1.16 


5.99 


95.12 


0.48 


9.02 


9.55 


Eridanus 


10.5 


4.26 


5.00 


1.03 


6.35 


77.01 


0.06 


9.32 


9.70 


Pal3 


11.0 


4.48 


3.40 


0.74 


12.74 


89.72 


0.04 


9.95 


9.95 


Pal4 


10.5 


4.62 


3.60 


0.77 


12.59 


93.77 


0.18 


9.96 


10.01 



iMcLaughlin & van der Marel (120051 we adopted an average 
value of 1.9 for Eridanus, which is not included in that work). 
The resulting values of f lc and f,/, are shorter than the age of 
the clusters (Table |2). Hence, some degree of central mass 
segregation is expected in the program clusters. This result is 
at odds with what found in the two outer halo GCs previously 
investigated, NGC 2419 and Pal 14, whose t rh (~ 18 and ~ 20 
Gyr, respectively) are lon ger than the clusters' age s (~ 12 and 
~ 10.5, respectively; see Dalessandro et al. 2009, and D08). 

4. THE BSS RADIAL DISTRIBUTION 

As de monstrated by previo us resul ts ( see e.g. iFerraro et all 
20041 iCarraro & Seleznevl l20TTt lSannaetal.1 120 121 
Salina setaLl 120121 and references therein), BSSs in 



most Galactic GCs are more centrally concentrated than 
normal stars. Since BSSs are more massive than the average, 
and since the half-mass relaxation time in those systems is 
much smaller than their age, this result is generally ascribed 
to the effect of dyn amical friction and mass segregation (e.g. 
Mape lli et alJl2006h . 

The clusters studied in this paper contain a significant num- 
ber of candidate BSSs (see Figured see also Sandquist 2005, 
hereafter S05). In order to investigate their radial distribution, 
the BSS populations have been selected in a homogeneous 
way in all clusters. We assumed the TO as the reference point 
and selected as bona-fide BSSs all the stars brighter than Vto 
and 3 x ay/ bluer than (V-I)to, where Vto ar >d {V-I)to 
are, respectively, the magnitude and colour of the TO as de- 
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FIG. 4. — Population selection in the four clusters. The bona-fide BSSs, 
located within the marked selection boxes, are shown as open circles. The 
reference population of RGB stars is marked with asterisks. 



TABLE 3 

Number of BSSs and RGB stars selected in each cluster. 



Cluster 


Nbss 


Nrgb 


AM 1 


31 


352 


Eridanus 


19 


187 


Pal 3 


15 


191 


Pal 4 


16 


354 



fined in Section 12.11 while ayj is the combined photometric 
uncertainty in the colour. The adopted selection boxes and 
the resulting bo na-fide BSS sample s are shown in Figure 2] 

According to lFerraro et al.l (11993b . in order to study the ra- 
dial distribution of BSSs, it is necessary to define a reference 
population which is expected to follow the cluster light distri- 
bution. Since the number of horizontal branch stars in these 
GCs is very low, we decided to adopt the RGB as the reference 
population. We selected as giants all the stars in the range of 
magnitude Vjo >V> Vro — 2.5, lying at a distance smaller 
than 3cr from the mean ridge line (see asterisks in Figure 2), 
where a is the uncertainty associated with the mean ridge line 
(horizontal bars in Figure [1}. This choice allows us to obtain 
a populous sample of reference stars in the same range of V 
magnitude, i.e., in the same condition of completeness as for 
the BSSs. The total number of BSSs and RGB stars in each 
cluster is listed in Table [3] 

4.1. Cumulative radial distribution and population ratios 

In Figure [5] we show the cumulative radial distribution of 
BSSs (solid lines) and RGB stars (dashed lines) as a func- 
tion of the projected distance from the cluster centre (r) nor- 
malized to the core radius r c estimated in Section 13.11 A 
Kolmogorov-Smirnov test has been performed to check the 
probability that BSSs and RGB stars are extracted from the 
same parent distribution. Based on the probability values (p) 
obtained for each cluster and reported in Figure |5] we con- 



FlG. 5. — Cumulative radial distributions of BBSs (solid lines) and RGB 
stars (dashed lines) as a function of their projected distance from the cluster 
centre normalised to the core radius (r/r c ), for each of the four clusters. The 
probability that the two populations are extracted from the same distribution 
is quoted in each panel. 

elude that BSS are more concentrated than the RGB in AM 
1 and Eridanus, while the result is inconclusive for Pal 3 and 
Pal 4. 

As second step to investigate the radial distribution of BSSs 
in the four clusters, we calculated the ratio between the num- 
ber of BSS (N B ss) and that of RGB stars (N RG b), counted 
in concentric annuli centred on the clusters' C grav . The radial 
distribution of Nbss /Nrgb in the four clusters is shown in Fig- 
ure|6] Consistently with what was found above, BSSs appear 
to be systematically more centrally concentrated than RGB 
stars. It is important to emphasise that, because of the insuf- 
ficient quality of the wide-field data, we are forced to limit 
our analysis to the region sampled by the HST observations, 
which do not reach the tidal radius of the program clusters. 
Thus, unfortunately we cannot characterise the full shape of 
the BSS radial distribution, and we are not able to conclude 
whether the distributions are bimodal or not. However this 
analysis demonstrates that dynamical friction was already ef- 
fective in centrally segregating BSSs in the program clusters. 

4.2. The Minimum Spanning Tree 

Given the relatively small number of BSSs, we applied 
the method of th e Minimum Spanning Tree (MST; see 
I Allison et al.l 120091 and references therein) as an alternative 
test to evaluate the degree of BSS segregation. The MST is 
the unique set of straight lines ("edges") connecting a given 
sample of points ("vertices"; in this case the star coordinates) 
without closed loops, such that the sum of the edge lengths 
is the minimum possible. Hence, the length of the MST 
is a m easure of the compactness of a given sample of ver- 
tices. ICartwright & Whitworthl (12004 showed that the de- 
gree of mass segregation in a star cluster can be measured 
by comparing the length of the MS T of two populations wit h 
different average masses (see also Schmeia & Klessen 20061) . 
Notice that the MST method is independent of the catalogue 
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TABLE 4 

Values of the MST tree (r M sT> and its la- uncertainty (Ar M sT) for the 

BSS IN THE FOUR TARGET CLUSTERS AND FOR NGC 2419 AND PAL 14. 



FIG. 6. — Relative frequency of BSS-to-RGB stars, as a function of the 
projected distance from the centre normalized to the core radius (r/r c ) for 
each of the four clusters. 

astrometric accuracy and of any assumption about the cluster 
center. Moreover, it gives quantitative measures of both the 
degree of mass segregation and its associated significance. 

Here we provide a detailed analysis of the segregation of 
BSSs with respect to RGB stars in the four program clusters, 
based on a recent modification of the MST method, where, in- 
stead of the di rect sum of n edges of length e, , their geometric 
mean is used (Olczak et al. 20 lib: 



7mst : 




: exp 



1 " 

- Vine,- 

n L — ^ 



(1) 



Following equation (2), we have computed 7^st f° r tne most 
massive population (i.e., the observed samples of n + 1 BSSs 
in each cluster). Then, we randomly extracted m sets of n+ 1 
objects from the reference population (the combination of the 
BSS and the RGB samples), and we computed their mean and 
standard deviation, TjJst an d ^7mst respectively. Finally, the 
level of BSS segregation with respect to giant stars and its 
associated uncertainty have been estimated as: 



T _ TmST AT 1 
1 MST - mn „ , MST - 



(2) 



Tmst 7mst 

The geometrical mean has the very useful property of effec- 
tively damping contributions from extreme edge lengths, so 
that Tmst of a compact configuration of even a few stars will 
not be much affected by an "outlier". Note that 7mst has the 
dimension of a length, while Tmst is dimensionless. A value 
of Tmst = 1 is found if the two populations have the same 
radial distribution, while Tmst > 1 indicates the presence of 
mass segregation meaning a more concentrated distribution of 
massive stars compared to the reference population. In the lat- 
ter case, the significance of mass segregation is provided by 
the value of k for which Tmst _ k ATmst = 1 ■ 

Table [4] lists the values of Tmst and the related la uncer- 
tainties computed for the BSS populations in the four tar- 



Cluster 


Tmst 


AFmsT 


AM 1 


1.71 


0.23 


Eridanus 


1.69 


0.27 


Pal 3 


1.56 


0.27 


Pal 4 


1.53 


0.29 


NGC 2419 


1.01 


0.07 


Pal 14 


1.07 


0.17 



get clusters, and, for comparison, for NGC 2419 and Pal 
14. In agreement with the results of Dalessa ndro et al.l (2008 ) 
and iBeccari et all (1201 lh . we find that BSS and RGB stars 
share the same radial distribution in NGC 2419 and Pal 14 
(Tmst — 1). Conversely, BSSs appear to be significantly more 
segregated than giant stars in the other four GCs. In particu- 
lar, AM 1 and Eridanus show the highest level of mass seg- 
regation (with Tmst — 1-7), with high statistical significance 
(2.5 — 3 a). Pal 3 and Pal 4 show a slightly lower degree of 
BSS segregation (Tmst — 1-5), at the ~ 2er level. 

5. SUMMARY 

Using HST-WFPC2 data we studied the BSS population in 
the central regions of four Galactic halo GCs, namely AM 
1, Eridanus, Pal 3 and Pal 4. These clusters, together with 
NGC 2419 and Pal 14, represent the entire group of GCs at 
R cc > 50 kpc. 

A proper comparison of the derived CMDs with theoretical 
isochrones confirms the younger ages (10.5-11 Gyr) of these 
systems with respect to the inner-halo GCs (see Figure [T). 
The impressive similarity of their RGB morphology (Fig. |2) 
also indicates a small metallicity difference (—1.41 <[Fe/H]< 
-1.58, depending on the cluster), in agreement with previous 
findings by S99 and D08. 

By complementing HST data with wide-field catalogues 
from ground based imaging, we have derived the most ex- 
tended radial density profiles from resolved star counts ever 
published for these clusters (Fig. |3J. All profiles are well fit 
by isotropic King models with the structural parameters listed 
in Table |2 

Unluckily, the insufficient quality of the wide field cata- 
logues did not allow us to study the BSS populations in the 
cluster outskirts and look for possible signatures of an exter- 
nal rising branch, similar to that found in m ost of the previ- 
ously surveyed GCs (see Sa nna et all 120121 and references 
therein). We therefore limited the analysis to the area covered 
by the HST data. BSSs have been selected in a homogenous 
way in the four program clusters and their radial distribution 
has been compared to that of RGB stars, taken as a reference. 
We found that BSSs are significantly more centrally concen- 
trated than giants in all four systems (see Figs. |5]|6]and Table 
ID). Since BSSs are assumed to be more massive than normal 
cluster stars, their higher central concentration is interpreted 
in terms of an effect of mass segregation. Indeed, the half- 
mass relaxation time of the four program GCs is found to be 
smaller than their age (see Table |2), thus indicating that dy- 
namical friction has been already effective in segregating the 
most massive stars in the cores. By measuring the degree of 
mass segregation with the Tmst test, we conclude that AM 1 
and Eridanus are the dynamically oldest systems in the group, 
while Pal 3 and Pal 4 are slightly less dynamically evolved, 
and NGC 2419 and Pal 14 show no signature of mass seg- 
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regati on yet (see also iDalessandro et al J [20081: iBeccari et alj 
1201 ll) . The results shown in this paper once more indicate 
that, indeed, BSSs represent the ideal population to investi- 
gate the dynamical state of dense stellar systems. 
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